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£H ' Abstract 

Aims. We investigate statistical and individual astrophysical properties of active galactic nuclei (AGNs), such as parsec-scale flux 
density, core dominance, angular and linear sizes, maximum observed brightness temperatures of VLBI core components, spectral 
index distributions for core and jet components, and evolution of brightness temperature along the jets. Furthermore, we statistically 
compare core flux densities and brightness temperature as well as jet spectral indices of y-ray bright and weak sources. 
Methods. We used 19 very long baseline interferometry (VLBI) observing sessions carried out simultaneously at 2.3 GHz and 8.6 GHz 
with the participation of 10 Very Long Baseline Array (VLBA) stations and up to 10 additional geodetic telescopes. The observations 
{J ; span the period 1998-2003. 

Results. We present here single-epoch results from high-resolution radio observations of 370 AGNs. Our VLBI images at 2.3 GHz 
and 8.6 GHz as well as Gaussian models are presented and analyzed. At least one-fourth of the cores are completely unresolved on 
the longest baselines of the global VLBI observations. The VLBI core components are partially opaque with the median value of 

Q ' spectral index of a core ~ 0.3, while the jet features are usually optically thin cr^, 0.7. The spectral index typically decreases along 

S_t , the jet ridge line owing to the spectral aging, with a median value of -0.05 mas -1 . Brightness temperatures are found to be affected by 

■ Doppler boosting and reach up to ~10 13 K with a median of ~2.5 x 10" K at both frequencies. The brightness temperature gradients 
along the jets typically follow a power law T b oc r~ 21 at both frequencies. We find that 147 sources (40%) positionally associated with 
y-ray detections from the Fermi LAT Second Source Catalog have higher core flux densities and brightness temperatures, and are 
characterized by the less steep radio spectrum of the optically thin jet emission. 

(N : 
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■ 1. Introduction from the NRAO archive^]. Since the project was initially focused 



in 



on performing geodesy and astrometry, the primary type of in- 



lO . m. i * urD i ■ * dfyw /x> u q i-. i * formation gathered was in the measured phases of the regis- 

. ' The long-term VLBI project RDV (Research & Development , . f , „ ,,,,, 

■ ,„ D n ■ , ■ j . e , .• f , ■ .* a . tered signals, whereas for us both the phases and amplitudes are 
, - VLBA) is designed to perform observations of bright, flat- „ • r ■ , , • , • 

v-N , . ,. T , equally important for restoring and analyzing the images. We 

_^ ■ spectrum, compact extragalactic radio sources. The project M , , „ r , , ° T _ . / & , • , , 

. f , ■ mm j tU j* r xt a c a a xTr>A/-v contacted the staff of the non- VLBA stations and received valu- 

• • started in 1997 under the coordination of NASA and NRAO. , , . . . _ T , T _. 

>' ™, ■ ,. , .. a oar^us able information about system equivalent flux densities (SEFD) 

, The simultaneous observations at 2.3 GHz and 8.6 GHz are car- , ,m r , • 

I j . ■ a ,li i- c ■ t u and system temperatures at different epochs of observations. We 

k> ' ned out bi-monthly, making up five to six sessions per year with \ . F . ,,,,,, 

K'S ■ tt, ^- • t ii » o/ „ TDA . a + + note that our imaging results should help to complete astrom- 

«_j the participation of all ten 25 m VLBA antennas and up to ten .. & , fe <.-.« ■. •• 

!3 ' j .• : .. f t , . rrr- . , |h nn nk tu • etry, i.e., improve the accuracy of the absolute positions of the 

■ geodetic stations (for detail see Petrov et al. 2009). The partici- ■" ,. . . . .„. , 

. , . D An , c \ e - N sources taking into account their milharcsecond structure, 

pation of the southern antennas such as HartRAO (South Africa) & 

and TIGO (Chile) allowed us to improve the Mv-coverage of low- In each experiment, about 100 active galactic nuclei are 

declination sources. One of the initial goals of the project was the scheduled. The total number of sources within RDV project is 

estimation of the precise absolute positions of compact sources currently approaching 1000. The main selection criteria for the 

and the impro vement o f the International Celestial Reference RDV targets x^re brightness and compactness of the sources 

Frame (ICRF; iMa et al.irT998l). w hich was extended by the in- d Petrov etal. | |2009|). In this paper, we present results of our 

elusion of 776 sources JFev et aTll2C)0l) an d later 3414 ra dio as- lma 8 ln g and analysis of raw archival RDV data for a sample 

tronomical objects constituting the ICRF2 dBoboltz et al.ll2010\ of 370 ob J ects ( Fl §- © observed during 19 sessions performed 

The second goal was to perform to precise geodesy and the de- Wlthln the P enod 1998-2003. If a source had been observed 

termination of antenna reference point change s associated with more than once ' we selected the e P och at whlch the dynamic 

the motion of tectonic platforms (iPetrov et aLll2009h . ran 8 e of the lma 8 e was lts hl g hesL The sam P le 18 dominated by 

quasars, with the weak-lined BL Lacs and radio galaxies making 

We have successfully used the same data for astrophysi- 



cal studies of active galactic nuclei (AGNs), drawing on data 1 http://archive.nrao.edu 
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up 8.3% and 7.8% of the sample, respectively. We used all 370 
sources available from the reduced observations for the analy- 
sis, because as shown by Lister & Marschen (1 19971 ). the most 
effective approach to studying Doppler-boosted sources such 
as AGNs is based on the investigation of samples with large 
number of objects. The selection effect may otherwise lead to 
incorrect conclusions about AGN properties as a class of ra- 
dio sources. Chronologically, the RDV project was one of the 
first among other large VLBI imaging surveys covering alto- 
gether the frequency range from 2 GHz to 86 GHz and in- 
cluding the VLBA Ca l ibrator Survey (VC S :lBeaslev et al.ll2002t 
Fomalont et all 120031; IPetrov et al.l 120051. l2006t iKovalev et alj 
20071: iPetrov et al.ll2008l), the VLBI Space Observato ry Program 



(VSOP: iFomalont et alj|205ot iDodson et all [20081). the VLBI 



Imaging and Polarime try Survey (VIPS: iHelmboldt et al.ll2007l 
IPetrov & Tavloril201 ll). the C altech-Jodrell Bank Flat-Spectrum 
sample (CJF; Pollack et al. 200 3J), t he co mplete sample of 
BL Lac ertae objects dGabuzda et al.1 1 2000|) . the VLBA 2 cm 
Survey (iKellermann et al.l 1998b . and the Monitoring Of Jets 
in Active galactic nuclei with VLBA Experiments (MOJAVE; 
iLister et al1l2009ah . the VLBI E xploration of Radio Astrom etry 
(VERA) and VLBA K-band dPetrov et alJ I2007L l2012h, the 



International Celestial R eference Frame (I CRF; Lanvi et all 
120 lOt ICharlot et all 1201 Oh and its extensions (iFev et al.ll2004l). 



and th e Coordinated Millimetre VLBI Array (CMVA: lLee et al.l 
[2008b surveys. 

First imaging results for early five 24-hour pre-RDV ses- 
sions with t he participation o f on ly the VLBA antennas wer e 
presented bv lFevetal.l (1 19961) and lFev & Charlotl (1 19971 l2000b . 
The im aging of some othe r RDV epochs prior to 1998 was also 
done bv lPiner et al.l(l2007l) . In this paper, we present and analyze 
our images from 19 experiments, the VLBI maps for the major- 
ity of RDV experiments incl uding maps presented in this paper 
are made available onlin dWl The astrome tric suitability of 80% 
of the ICRF sources was investigated by ICharlot et al.l (120081) . 
The RDV observations proved to be successful for spectral in- 
dex studies including the opacity effect detec tions reported by 
IKovalev et alj d2008l) and later confirmed by ISokolovsky et alj 
( 20111) . Another application is the measurement of the multi- 
epoch 8.6 GHz apparent speeds of the jet components. The ini- 
tial kinematics results for 54 sources based on 8.6 GHz VLBI 
images in Radio Reference Frame Image Databas e from the 
first fiv e years (1994-1998) of data were obtained bv lPiner et alJ 
d2007l) . Extending the period of observations up to ten years and 
enlarging a number of epochs up to 50 will both increase the 
accuracy of the apparent speed me asurements and a llow us to 
investigate possible jet acceleration dPiner et al 1 120121) . 

The observational material of the RDV project has a number 
of advantages compared with the data from many large VLBA 
experiments cited above: (i) improved coverage of Mv-plane; (ii) 
simultaneous dual-frequency mode; (iii) high time-sampling fre- 
quency (nearly every two months) for a subsample of about 60 
sources. This thir d feature permits the study of jet kinematics 
dPiner et al ]l2007h . while the other two are important for the pur- 
poses of this paper, in which we discuss the compactness of the 
sources and both their spectral indices and brightness tempera- 
tures, and compare these properties for y-ray bright and weak 
sources. 

Among the 370 RDV sources presented in this paper, 
147 objects have been positionally cross-identified with the 
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2 http://astrogeo.org/images 

3 http://rorf . usno.navy.mil/rrf id. shtml 

4 http : //www . obs . u-bordeauxl . fr/BVID 



Right Ascension 

Figure 1. Sky distribution of 370 observed sources in Aitoff 
equal-area projection of the celestial sphere in equatorial coor- 
dinates. The thick line is the Galactic equator. 



high-confidence (>4<x) y-ray detections associated with known 
AGNs, according to the Fermi Large Are a Telescope Second 
Source Catalog (2FGL; lAbdo et aiTl201 ll) . These objects are 
both targets of follow-up observational VLBI sessions and 
prospects for establishing connections between y-ray emission 
and parsec-scale jets observed in the radio domain. 

The structure of this paper is as follows: in §[2] we describe 
the observations and data reduction technique, including the cal- 
ibrating, imaging, and model-fitting procedures; in §[3] we dis- 
cuss our results; and our main conclusions are summarized in 
§[4] Throughout the paper, the ACDM cosmological model with 
H = 70 km &- 1 Mpc-\ Q„, = 0.3, and £l A = 0.7 is adopted. The 
spectral index a is defined in the convention S oc y". All position 
angles are given in degrees east of north. We use the term "core" 
as the apparent origin of AGN jets that commonl y appears as the 
brightest feature in VLBI images of blazars (e.g jLobanovIl 99*81 ; 
lMarschej|2008l) . 



2. Source sample, observational data, and data 
processing 

2. 1 . Source sample and its completeness 

The sample of 370 sources includes 251 quasars (67.8%), 46 BL 
Lacertae objects (12.4%), 31 radio galaxies (8.4%), and 42 op- 
tically unidentified sources (11.4%). The redshifts are currently 
known for 306 objects (83%). The redshift distribution (Fig.[2]i 
ranges from Zmm = 0.00436 for the galaxy J 1230+ 1223 (M87) to 
^max = 4.715 for the quasar J1430+4204 with a median value of 
the distribution close to z = 1 . The general characteristics of the 
sources such as object name, its alternative name, coordinates of 
J2000.0 right ascension and declination, optical class, redshift, 
and linear scale in parsecs per mas are listed in Table [T] 

To investigate the completeness of the studied sample, we 
determined a cumulative log - log S dependence (as shown 
in Fig. [3]) for the studied RDV subsample and compared it with 
the Radio Fundamental Catalog (RFC) in the matched sky area, 
with declinations above -45°. The RFC provides (i) precise po- 
sitions with milli-arcsecond accuracies, (ii) estimates of corre- 
lated flux densities at baselines from 1000 km to 8000 km, and 
(iii) maps for thousands of compact radio sources produced by 
the analysis of all available VLBI observations under absolute 
astrometry and geodesy programs. The RFC sample constructed 
from a blind survey is statistically flux density complete down to 



|http:// astrogeo . org /rfc| 
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Figure 2. Redshift distribution of 306 sources with the median 
and maximum values of 0.98 and 4.715, respectively. 



200 mJy dKovalev et alj|2007h IPetrov et al.ll2008l) . We note that 
the difference of 1-2 sources seen for the brightest objects is 
due to the flux density variability and that the RFC provides the 
flux densities averaged over all available epochs, while the RDV 
statistics is based on single-epoch measurements. As seen from 
Fig. |3] the studied RDV subsample is complete down to ~ 1.5 Jy. 
At lower flux densities, the sample is statistically incomplete but 
representative of compact A GNs owing to the o riginal selection 
strategy of the RDV project (IPetrov et al.ll2009l) . 

2.2. Observational data 

We reduced and present here single-epoch imaging results of 19 
observational 24-hour sessions (RDV11, 13, 15, 18, 21, 24, 27, 
29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, and 41) carried out 
in a period between October 1998 and September 2003 using a 
global VLBI array of 18 to 20 stations including 10 VLB A an- 
tennas and 14 geodetic and EVN radio telescopes (Tabled ca- 
pable of recording the VLBA modes, using 8 to 10 non-VLBA 
stations in each session. The dates of observations are listed in 
the first column of Table |4] The position s for all ante n nas gi ven 
in geocentric coordinates are listed in IPetrov et al.l (120091) of 
Table 3. The observations were performed in a standard dual- 
frequency geodetic mode, registering the signal in right circular 
polarization simultaneously at 2.3 GHz (S-band) and 8.6 GHz 
(X-band), and using 1-bit sampling. Each band was separated 
into four 8 MHz (4 MHz for RDV11) intermediate frequencies 
(IFs) giving a total observing bandwidth of 32 MHz (16 MHz for 
RDV1 1) at each band. The IFs spanned 140 MHz at S-band and 
490 MHz at X-band (Table©. The data were correlated with the 
VLBA correlator at the Array Operation Center in Socorro with 
4 s accumulation periods. Most sources were observed during six 
scans of six minutes each, corresponding to a total tracking time 
of about 40 min. A few sources were observed in up to 15 scans, 
making to a total on-source time of about 60 min. The scans are 
scheduled over a number of different hour angles to maximize 
the (u, v) plane coverage. The median number of independent 
measurements of visibility function on individual interferometer 
baselines averaged within a scan is about 200. No special cali- 
brators were scheduled since all the target sources were compact 
and bright. More than 92% of the sources in 19 processed ses- 
sions have VLBI fluxes greater than 200 mJy (see Fig. [9] top 
panels). 
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Figure 3. Cumulative log N - log S dependence at 8.6 GHz for 
the studied RDV sample of 370 sources (open circles) in com- 
parison with the Radio Fundamental Catalog (stars). The dotted 
line indicates a completeness level of about 1 .5 Jy for the studied 
sample. 



Observations were scheduled usi ng the automa tic mode of 
the VLBI scheduling software sked (lGipsonll20ldh with some 
subsequent improvements to maximize mutual visibility, ensure 
adequate sky coverage, minimize slewing time, etc. The detailed 
description of the sc heduling of the glo bal geodetic VLBI ses- 
sions is discussed bv lPetrov et all d2009l) . 

The special feature of these observations is the use of sub- 
arrays. The involvement of the large number of stations located 
on different continents prevents the antennas from being able to 
simultaneously track the same target source altogether. To max- 
imize the on-source time at each station during the whole ex- 
periment, the array of radio telescopes is divided into several 
subarrays to ensure that a source is visible for every station in 
a corresponding subarray during a particular scan. The number 
of subarrays ranges from 6 to 8 in the reduced experiments. We 
note that each antenna is included not only in some particular 
subarray, but in several ones for different scans. Thus, the major- 
ity of sources are observed with almost all antennas participating 
in a session. The exceptions are the high-declination sources for 
the southern stations and vice versa. 

The 5 m VLBI antenna MV-3 GGAO near Washington also 
took part in the observations but we did not use the data from 
this site because of its extremely high noise level (SEFD > 
5 x 10 4 Jy). In addition, GGAO in a number of experiments 
was tagged along after the schedule was generated using other 
stations in a session. GGAO was added to each scan that the an- 
tenna could slew to in time, and the scan length for GGAO was 
set to the maximum of all the stations already participating in the 
scan. This resulted in a significant number of non-detections on 
the GGAO baselines. 

The southern hemisphere antennas, the 26 m radio telescope 
HartRAO in South Africa (latitude of -26°) and the 6 m TIGO 
antenna in Concepcion, Chile (latitude of -36°) have partici- 
pated in 15 out of the 19 analyzed observational sessions. This 
significantly increased the resolution in the north-south direction 
and allowed us to image low-declination sources with almost cir- 
cular restoring beam. The contribution of data from even such a 
small antenna as TIGO can significantly improve the quality of 
the final image (Fig.|4|, if the source is bright (S > 0.5 Jy). 

Increasing the number of antennas in the array from 10 
(VLBA) to 18-20 (global VLBI) allowed us on-average to: 
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Table 1. Parameters of observed sources. 



Source 


IVS name 


Right ascension 


Declination 


Optical 


z 


Scale 


2FGL 






J2000.0 




J2000.0 


class 




pc mas 


member 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


J0006-0623 


0003-066 


00 06 13.892887 


-06 


23 35.33545 


B 


0.347 


4.913 




JUU1 1— ZOlZ 


UUUo— 2 04 


UO It U1.Z40/3/ 


-26 


12 33.37719 





1 ooq 
I.W5 


o.lu3 




J0017+8135 


0014+813 


00 17 08.474901 


+81 


35 08.13655 


Q 


3.387 


7.407 




J0019+7327 


0016+731 


00 19 45.786379 


+73 


27 30.01753 


Q 


1.781 


8.450 




J0022+0608 


0019+058 


00 22 32.441214 


+06 


08 04.26889 


B 






Y 


J0027+5958 


0024+597 


00 27 03.286191 


+59 


58 52.95918 


U 








J0035+6130 


0032+612 


00 35 25.310644 


+61 


30 30.76129 


u 








J0050-0929 


0048-097 


00 50 41.317382 


-09 


29 05.21043 


B 






Y 


J0059+0006 


0056-001 


00 59 05.514929 


+00 


06 51.62077 


Q 


0.717 


7.215 




J0102+5824 


0059+581 


01 02 45.762379 


+58 


24 11.13659 


Q 


0.644" 


6.902 


Y 



Columns are as follows: (1) J2000.0 IAU name; (2) IVS name (B 1950.0 IAU name or alias); (3) right asc ension (12000.0) in hours, mi nutes, and 
seconds; (4) declination (J2000.0) in degrees, minutes, and seconds; (5) optical classification according t o | Verpn-Cettv & Ve ron (2006), where Q 
is a quasar, B is a BL Lace rtae object, G is an active galax y, and U is unidentified, (6) red shift as gi ven bv |Veron-C ettv & Ve ron ( 2006 ), the a flag 
indica tes the red shift fromlSowards-Emmerd et all (120051), b - fromlCarilli et al.l dl998|), c - fromlDrake et all ( 120041). d - from | Afanasiev et al.l 
( 12003 1 ). e - fromlvon Montignv et al.l dl995|). f - fromlNilsson et all (l2008h . g - from lBest et all d2003h . h - from lHewitt & Burbided dl989h . k - 
from lSnellen et aL I d2002bl) , and / - from Sne llen et al.l d2002al); (7) an gular scaling conversion in parsecs per milliarcsecond; (8) membership in 
the Fermi LAT Second Source Catalog according to lAbdo et alj <|20 1 lh . The coordinate positions are taken from the Radio F undamental Catalog 
( http ://astrogeo . org/rf c ( as derived from an analysis of all VLBI ob s ervations made in absolute astrometry and geodesy mode (Beasley et al. 2002; 
Fomalont et al. 2003; Petrov et al. 2005, 2006: Kovalev et al. 2007; Petrov et al. 2008). Table[T]is published in its entirety in the electronic version 
of the Astronomy & Astrophysics. A portion is shown here for guidance regarding its form and content. 



Table 2. Antennas participating in the RDV observations. 



Station 


Location 


Code Diameter 






m 



46 
25 
25 
26 
20 
25 
26 
25 
20 
25 
20 
25 
32 
25 
32 
20 
20 
25 
25 
25 
6 
32 
18 
20 



- increase the number of possible baselines by a factor of ~2.5, 

- increase the angular resolution by a factor of ~1.5, 

- decrease the noise level of the resultant VLBI image by a 
factor of ~2. 

The maximum projected baseline of about 12 300 km between 
Owens Valley (VLBA) and HartRAO that was realized in the 
RDV experiments is close to the limit for Earth-based VLBI ob- 
servations. In Fig. [5] we compare the typical (u, v) coverage at 
8.6 GHz provided by the VLBA and the whole array used for 



Table 3. The range of used frequencies. 



IF Frequency range^, [MHz] 



1 


2232.99 


- 2240.99 


2 


2262.99 


- 2270.99 


3 


2352.99 


- 2360.99 


4 


2372.99 


- 2380.99 


5 


8405.99 


- 8413.99 


6 


8475.99 


- 8483.99 


7 


8790.99 


- 8798.99 


8 


8895.99 


- 8903.99 



" S-band IFs were shifted up 12 MHz since the RDV27 session 
(2001.04.09) to get further away from direct broadcast satellite band. 



high-, medium-, and low-declination sources. The correspond- 
ing spatial frequency plane coverages at 2.3 GHz are virtually 
identical but scaled accordingly. 

2.3. Data processing 

Initial calibration was performed w ith the NRAO Astronomical 
Image Processing System (AIPS: iGreisenl 120031) by applying 
techniques adopted for subarrayed data sets. The individual IFs 
were processed separately throughout the data reduction, and 
were averaged together when making the final images. The an- 
tenna gains and system temperatures measured during the ses- 
sions were used for the amplitude calibration. For some of the 
non-VLBA stations, there was no direct information about the 
antenna's gain at the epoch of an experiment. In these cases, we 
estimated the degrees per flux unit (DPFU) parameter (gain) for 
each band using the information provided by the staff of the ob- 
servatories about the SEFD measured on a date closest to the 
epoch of observations (to within several days) and the system 
temperature averaged over the current session. Global gain cor- 
rection factors for each station for each IF were derived from 
the statistical results of self-calibration for all the sources within 



Algonquin Park 


Canada 


AP 


Brewster 


USA 


BR 


Fort Davis 


USA 


FD 


Gilcreek 


USA 


GC 


Green Bank 


USA 


GN 


Hancock 


USA 


HN 


HartRAO 


South Africa 


HH 


Kitt Peak 


USA 


KP 


Kokkee 


USA 


KK 


Los Alamos 


USA 


LA 


Matera 


Italy 


MA 


Mauna Kea 


USA 


MK 


Medicina 


Italy 


MC 


North Liberty 


USA 


NL 


Noto 


Italy 


NT 


Ny Alesund 


Norway 


NY 


Onsala 


Sweden 


ON 


Owens Valley 


USA 


OV 


Pie Town 


USA 


PT 


St. Croix 


USA 


SC 


Concepcion 


Chile 


TC 


Tsukuba 


Japan 


TS 


Westford 


USA 


WF 


Wettzell 


Germany 


WZ 
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J0017+8135 at 8.6 GHz in RR 2002 Mar 06 




/.--.'i '• 9- iff 



J1310+3220 at 8.6 GHz in RR 2002 Jan 16 



J0457-2324 at 8.6 GHz in RR 2003 Mar 12 



<iir 



200 -200 200 -200 200 -200 

U (10 6 A) U (10 6 A) U (10 6 A) 

Figure5. Coverages of the («, v) plane at 8.6 GHz for a high-declination (J0017+8135, left), medium-declination (J1310+3220, 
middle), and low-declination (J0457-2324, right) source with the baseline projections produced by the VLBA antennas only (red 
dots) and the global VLBI array (red + blue dots). 



Table 4. Observing log with non-VLBA telescope gains. VLBA antennas participated in all sessions. 



Date 2.3/8.6 GHz DPFU", [mK Jy~'] 





AP 


GC 


GN 


HH 


KK 


MA 


MC 


NY 


NT 


ON 


TC 


TS 


WF 


WZ 


1998/10/01 




85/125 


78/70 




70/41 




140/128 


24/69 




23/45 






43/39 




1999/03/08 




95/67 


49/76 


33/21 


60/47 




138/152 


30/77 




24/48 






47/40 


63/68 


1999/05/10 




119/74 




36/38 


60/43 




131/143 


31/72 




23/34 






43/36 


66/79 


1999/12/20 




69/84 


50/72 


38/30 


60/47 




140/146 


30/77 




23/45 




176/269 


45/38 


57/38 


2000/05/22 




88/74 




36/24 


67/39 


27/45 


128/118 


34/73 








223/271 


46/37 


71/80 


2000/12/04 




82/94 




36/22 


65/36 


47/49 


180/134 


32/76 








205/269 


44/39 


68/82 


2001/04/09 




95/80 
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129/83 
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309/122 
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2002/12/11 
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59/40 


32/41 
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28/34 


5/8 
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53/44 


29/45 
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4/8 
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61/36 
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2003/06/18 










57/42 


-/47 








31/39 


6/9 


220/183 


66/31 


-/78 


2003/09/17 




86/74 






61/46 






26/54 




20/36 




175/171 


27/37 


75/95 



" Efficiency (degrees per flux unit) of a radio telescope in right circular polarization in zenith direction. In all cases, the gain curve is assumed to 
be flat. 



an experiment taking into account the well-measured VLBA an- 
tenna gains. If the median value of the gain correction distribu- 
tion deviated by more that 10% from 1, the whole calibration 
procedure was performed again, where we applied the corre- 
sponding corrections by running the AIPS task CLCOR prior 
to the stage of phase calibration. We estimated the accuracy of 
amplitude calibration to be at the level of ~10%. The non-VLBA 
station sensitivities obtained from the log-files and then defined 
more accurately from the gain correction technique are listed in 
Table |U 

Because the phase calibration tones were available for the 
VLBA stations only (non-VLBA stations have no phase-cal de- 
tectors), the phase corrections for station-based residual delays 
and delay rates were found and applied using the AIPS task 
FRING in two steps. First, the manual fringe fitting was run on 
a short interval of either one or two minutes on a bright compact 



source to determine the relative instrumental phase and resid- 
ual group delay for each individual IF. As the array of radio 
telescopes is divided into several subarrays, we had to run the 
procedure on a few different scans to find the solution for each 
antenna. Secondly, the global fringe fitting was run by specify- 
ing a 4 min solution interval, a point-like source model, and a 
signal-to-noise ratio (S/N) cutoff of 5 to omit noisy solutions. 
The fraction of inadequate solutions was typically at a level of 
(8-10)% at X-band and (10-13)% at S-band. 

Amplitude bandpass calibration was not applied because it 
had no significant effect on the dynamic range of our final im- 
ages. We used the AIPS task SPLIT to write the fully cali- 
brated data for each individual source to separate X- and S-band 
FITS-files averaging data in frequency within each of four IFs at 
each band, and making no averaging in time. First, phase self- 
calibration with a point-source model was done using the AIPS 
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Figure 4. Naturally weighted CLEAN images at 2.3 GHz of 
the low-declination source J0403-3605 without (left) and with 
the data (right) from TIGO antenna. The contours are shown at 
the same levels, starting from 2 mJy beam -1 and increasing by a 
power of 2. The axes of each image are given in milliarcseconds. 
The shaded ellipse in the lower left corner of each image indi- 
cates the full width at half maximum (FWHM) of the restoring 
beam. 



task CALIB, which provides: (i) the possibility of finding the 
first solution using all IFs, (ii) automatic flagging the data with 

low S/N ratio. 

19741). phase an d amplitude self- 



CLEA Ning dHogbom 



calibration (Jenniso 



ogbom 1^/41). phase an d amplitude selt- 
nll 1 958 :lTwiss et al.ll 19601). and hybrid imag- 



ing dReadhead et alJl980HSchwabll980HCornwell & Wilkinson! 
11981b were performed in the Caltech DIFMAP (IShepherdl 199% 
package using an au tomated approach suggested by G. Taylor 
dPearson et al.ll 19941) . which we adapted to the RDV data set. In 
all cases, a point-source model was used as an initial model for 
the iterative procedure. Final maps were produced by applying a 
natural weighting of the visibility function. The spanned band- 
width of the four IFs in each band is relatively small (6% of 
fractional bandwidth in both bands), thus no spectral correction 
technique was applied. In the model fitting also performed in 
DIFMAP, we used a minimum number of circular (and in some 
cases elliptical) Gaussian components that after being convolved 
with the restoring beam, adequately reproduce the constructed 
brightness distribution. 



3. Results 

3.1. Calibrated VLBI images 

Final, naturally weighted VLBI images at 2.3 and 8.6 GHz for 
370 sources are shown in Fig. [6] The dynamic range of the im- 
ages determined as a ratio of the peak flux density to the rms 
noise level ranges from 66 to 7042 with a median of ~1200 at 
X-band, and ranges from 106 to 4789 with a median of ~1000 
at S-band. The typical rms noise level is ~0.4 mJy beam -1 at 
X-band and ~0.5 mJy beam -1 at S-band. 

In Table |5j we summarize the VLBI map parameters: (1) 
source name in J2000.0 notation; for X-band: (2) peak flux den- 
sity in Jy beam -1 ; (3) lowest contour level in mJy beam -1 ; (4) 
flux density integrated over entire map in Jy; (5) flux density of 



the most compact fraction in Jy as measured from the longest 
baselines; (6) major axis of the full width at half maximum 
(FWHM) of the restoring beam in milliarcsec; (7) minor axis 
of FWHM of restoring beam in milliarcsec; (8) position angle 
of major axis of restoring beam in degrees; columns (9)— (15) 
are the same as (2)-(8) but for S-band. Flux densities shown in 
columns (4) and (11) were calculated as a sum of all CLEAN 
components over VLBI image. Columns (5) and (12) give the 
flux densities of the unresolved fractions 5 unres and are estimated 
as the median of the visibility function amplitudes measured at 
spatial (u, v) radius r uv = (u 2 + v 2 ) 1/2 > 270 MA for X-band and 
r uv > 72 MA for S-band, corresponding to the range between 0.8 
and 1.0 of the maximum projected spacing. For 15 sources at 
X-band and 6 sources at S-band, no estimates of unresolved flux 
density are given because no data were available on the baselines 
used for the calculations. 

In addition to Table [5] all calibrated data including flexi- 
ble image transport system (FITS) files of naturally weighted 
CLEAN images, FITS files with (u, v) data sets, images in 
postscript and GIF-formats, plots of correlated flux density ver- 
sus baseline length projection, and values of total VLBI flux den- 
sity for all the sources at all available epochs listed in Table|4]are 
accessible online through the HTML interface^ 

3.2. Parsec-scale structure modeling 

The observed sources have very complex three-dimensional 
structures on parsec and dekaparsec scales. However, since they 
are located at cosmological distances they provide us with the 
possibility of investigating only their two-dimensional bright- 
ness distributions in the plane of the sky, which can be mod- 
eled by a limited number of Gaussian components. The number 
of components is usually determined by reconciling two oppos- 
ing conditions: (i) it should be minimal; (ii) all the components 
are convolved with the restoring beam, adequately reproducing 
the general structure of the obtained brightness distribution. The 
typical morphology of the sources is represented by a one-sided 
core-jet structure, which is the result of str ong selection effect s 
and the Doppler boosting of jet emission (ICohen et al.ll2007h . 
Only 11 sources (3% of the sample) have two-sided structures 
detected in the S/X images. There are also extreme cases when 
(i) a source is highly compact without any indications of a jet, 
presumably because of lack of spatial resolution and/or sensitiv- 
ity of the interferometer at the observed frequency; (ii) a source 
has a prominent jet that is resolved even in a transverse direction. 
In Fig.EI we plot the images of J1935+2031 and J0738+1742 as 
examples. 

Structure modeling of all sources was performed with the 
procedure modelfit in DIFMAP package by fitting several circu- 
lar Gaussian components to the self-calibrated visibility data and 
minimizing^ 2 in the spatial frequency plane. Elliptical Gaussian 
components were adopted when the use of circular Gaussian 
component resulted in the appearance of axisymmetric noise on 
the residual map near the modeled component. The model fits 
listed in Table [6] provide flux densities, positions, and sizes for 
the bright and distinct jet features. All the positions are given 
with respect to the core component. The quality of model fit- 
ting is generally quite good as can be seen in Fig. [8] where we 
plot the total flux density of the Gaussian components of the fit 
and CLEANed flux density calculated as a sum of all CLEAN 
components from the image. 
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Figure 6. One-epoch 2.3 GHz and 8.6 GHz VLBI images of 370 active galactic nuclei. Image parameters are listed in Table [5] 
Gaussian models fitted to the visibility data at each frequency are listed in Table [6] The scale of each image is in milliarcseconds. 
The shaded ellipse in the lower left corner indicates the FWHM of the restoring beam (with natural weighting). The contours are 
plotted at increasing factors of 2 times the lowest contour drawn at 3<x level. A single negative contour equal to the bottom positive 
contour is also shown. A cell size of 0.2 and 0.5 mas per pixel was used for 8.6 and 2.3 GHz maps, respectively. Fig. |6]is available 
in its entirety in the electronic version of the Astronomy & Astrophysics. 
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Table 5. Parameters of VLBI maps. 



Source 


2.3 GHz 


8.6 GHz 




peak 


clev 


S VLBI 


lJ unres 








peak 


clev 


S VLBI 


lJ unres 




^maj 






Jy tar 1 


mjy tar 1 


Jy 


Jy 


mas 


mas 


deg 


Jy tar 1 


mJy tar 1 


Jy 


Jy 


mas 


mas 


deg 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ii) 


(12) 


(13) 


(14) 


(15) 


JUUUo— Uozj 


1.33 


3.12 


2.33 


0.66 


3.42 


8.03 


4 


1.41 


2.92 


2.07 


1.18 


0.91 


2.23 


4 


JOOl 1-2612 


0.39 


1.90 


0.41 


0.36 


2.06 


5.56 


-5 


0.31 


1.60 


0.34 


0.27 


0.56 


2.03 


-10 


J0017+8135 


0.52 


0.96 


0.76 


0.40 


2.55 


2.71 


54 


0.49 


0.76 


0.75 


0.29 


0.68 


0.69 


7 


JOOl 9+7327 


0.57 


1.20 


0.76 


0.35 


2.51 


2.86 


20 


0.15 


1.00 


0.34 


0.13 


0.68 


0.74 


19 


J0022+0608 


0.42 


0.79 


0.45 


0.40 


3.12 


6.50 


-2 


0.41 


0.58 


0.44 


0.35 


0.84 


1.76 


-1 


J0027+5958 


0.07 


0.90 


0.08 


0.07 


2.81 


3.63 


-23 


0.06 


0.90 


0.07 




0.65 


1.02 


-24 


J0035+6130 


0.07 


2.06 


0.11 


0.09 


4.01 


5.00 


19 


0.08 


1.94 


0.11 


0.07 


1.07 


1.34 


2 


J0050-0929 


0.53 


1.93 


0.58 


0.45 


2.97 


9.62 


-8 


0.70 


1.65 


0.75 


0.61 


0.80 


2.62 


-7 


J0059+0006 


0.42 


5.80 


1.17 


0.25 


4.38 


7.68 


-4 


0.16 


5.40 


0.43 


0.06 


1.14 


2.18 


-8 


J0102+5824 


0.61 


0.92 


0.81 


0.42 


2.63 


3.25 


32 


1.55 


0.88 


1.68 


1.46 


0.67 


0.87 


35 



Columns are as follows: (1) source name (J2000.0), (2) peak flux density in image, (3) lowest contour in image, (4) total flux density from VLBI 
map, (5) flux density from the most compact component on map, (6) FWHM minor axis of restoring beam, (7) FWHM major axis of restoring 
beam, (8) position angle of major axis of restoring beam, (9)— (15) the same as (2)-(8), respectively. Table [5] is published in its entirety in the 
electronic version of the Astronomy & Astrophysics. A portion is shown here for guidance regarding its form and content. 



Table 6. Source models. 



Source 


B 


s 


r 


<P 


Maj. 


Ratio 


P.A. 






mjy 


mas 


deg 


mas 




deg 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


J0006-0623 


S 


1527 
46 


2.96 


-55 


1.87 


0.93 
1.00 


-37 






772 


5.56 


-77 


2.79 


1.00 






X 


1551 






0.81 


0.34 


7 






112 


1.26 


-68 


0.50 


1.00 








115 


2.23 


-65 


1.58 


1.00 








128 


5.34 


-76 


2.77 


1.00 








180 


6.58 


-78 


0.87 


1.00 




JOOl 1-2612 


s 


401 






0.57 


1.00 






X 


327 






0.26 


1.00 





Columns are as follows: (1) source name (J2000.0), (2) frequency band, 
(3) fitted Gaussian flux density, (4) position offset from the core com- 
ponent, (5) position angle of the component with respect to the core 
component, (6) FWHM major axis of the fitted Gaussian, (7) axial ra- 
tio of the fitted Gaussian, (8) major axis position angle of the fitted 
Gaussian. Table [6] is published in its entirety in the electronic version 
of the Astronomy & Astrophysics. A portion is shown here for guidance 
regarding its form and content. 



3.3. Compactness 

One of the main goals of the RDV program was to construct 
the radio reference frame. The sources included in the observing 
program thus have to lie at cosmological distances, be bright and 
compact, i.e., to radiate the emission from the extremely limited 
sky area. The most suitable objects for achieving these goals are 
active galactic nuclei, which easily satisfy the aforementioned 
criteria. The compactness of blazars on arcsecond scales derived 
as SvLBi/Stot, where S tot is the single-dish flux densities and 
S vlbi is the integrated flux density from a VLBI image, is at a 
ve ry high level of ~0 . 9 at 2 GHz and 8 GHz as initially reported 
bv lPopov & Kovalev 1 19991) for a s mall s ample of 20 sources and 
later confirmed by Kovalev et al.l (120051) for a complete sample 
of 250 flat-spectrum AGNs observed with the VLBA at 15 GHz. 
This implies that nearly all of the emission from these sources 
is generated on the milliarcsecond scales, probed by VLBI ob- 



servations, and that the contribution from extended emission on 
kiloparsec scales is typically negligible. 

The compactness of the sources from our sample on mil- 
liarcsecond scales, determined as a ratio of S unres to S vlbi, dis- 
tributions of which are presented in Fig. |9] is statistically lower 
(Fig. [TUJ top panels) with median values of 0.51 at 2.3 and 
8.6 GHz (Table©. According to both K-S and T-test, the distri- 
butions are indistinguishable. At higher frequencies, the VLBI 
compactness is also at the similar level. Thus, 68% of the sam- 
ple of 250 extragalactic radio so urces observed at 15 GHz have 
compactnesses greater than 0.5 (iKovalev et al .1120051) . The me- 
dian value of co mpactness at 86 GHz is about 0.5 as reported by 
iLee etalJd2008l) . 

We note that the overwhelming majority of the sources are 
also strongly core dominated (Fig. [10] bottom panels). In 50% 
of sources, the emission from the core region accounts for more 
than 75% of the total VLBI flux. The information on S vlbi, 
S unres, and compactness will be useful for selecting the most suit- 
able candidates (highly-compact bright sour ces) to be studie d 
with the RadioAstron space VLBI mission (lKardashevlll997l) . 
At the same time, there is a fraction of the sources that manifest 
strong jets (M87, 3C 120, etc.), which are of particular interest 
for detailed astrophysical studies of their extended outflows on 
parsec scales. 

3.4. VLBI core properties 

We define VLBI core as a compact, bright emitting region at 
the narrow end of the jet, where it becomes optically thick. 
Typically, VLBI cores are often unresolved along one or both 
axes of the restoring beam and have a flat radio spectrum. Since 
the core is usually the brightest feature of the source it is placed 
by the phase self-calibration procedure at the center of the map. 
In such a case, the core component always has a very small but 
nonzero shift from the phase center because of the close pres- 
ence of other bright components down the jet. Significant offsets 
of the VLBI core components from the occur phase center of an 
image are rare and happen when one of the jet components is 
brighter than the core. This kind of structure has been found for 
18 sources at X-band, and for 27 ones at S-band. 
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Figure 7. Naturally weighted CLEAN images of J 1935 +2031 
at 8.6 GHz (top, dynamic range of 570) and J0738+1742 at 
2.3 GHz (bottom, dynamic range of 1750) as examples of highly 
compact and jet-resolved sources, respectively. The axes of each 
image are given in milliarcseconds. The shaded ellipse in the 
lower left corner indicates the FWHM of the restoring beam. 
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Figure 8. Total flux density from Gaussian fits versus total flux 
density from VLBI CLEAN image for the data at 2.3 GHz (left) 
and 8.6 GHz (right). 
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Figure 9. Distributions of the VLBI flux density S vlbi (top pan- 
els) and the flux density of the most compact component S umes 
(bottom panels). 
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Figure 10. Distributions of the compactness 5 unres /5vLBi (top 
panels) and the VLBI core dominance S CO ie/S vlbi (bottom pan- 
els). 



Model fitting enabled us to measure flux densities and sizes 
of the core components. The distributions of the core flux densi- 
ties are plotted in Fig. [TT] (top panels). Owing to the high com- 
pactness of the cores, they are often unresolved. We considered 
the core to be unresolved if the size of a circular Gaussian com- 



10 



A. B. Pushkarev, Y. Y. Kovalev: Single-epoch VLBI imaging of active galactic nuclei at 2 GHz and 8 GHz 



60 - 



40 



30 - 



2.3 GHz - - 



..6 GHz - 



60 



40 



20 



0.1 



10 0.1 

Score ( J y) 



60 - 



40 



20 - 



I Measured 
□ Limits _ - 



10 



I Measured 
D Limits 



80 



60 



40 



20 



1 



0.5 
, (mas) 



60 



40 



20 



I Measured 
D Limits 



1.5 



I Measured 
D Limits 



60 



40 



20 



10 20 30 



10 



Figure 11. Distributions of flux (top), angular size (middle), 
and projected linear size (bottom) for the core components at 
2.3 GHz (left) and 8.6 GHz (right). The unshaded boxes repre- 
sent upper limits for those sources with unresolved core compo- 
nent. For the core components fitted with elliptical Gaussians, 
the major axis is used for the angular core size. 
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Figure 12. Brightness temperature distribution of VLBI core 
components in source rest frame at 2.3 GHz (left) and 8.6 GHz 
(right) with median value of 2.51 x 10 11 K and 2.60 x 10 11 K, re- 
spectively. The unshaded boxes represent lower limits for those 
sources with either unknown redshift or unresolved core compo- 
nent. 



ponent or at least one of the axes of an elliptical Gaussian com- 
ponent was smaller than the respective resolution limit (which is 
a function of the restoring bea m parameters and S/N at the core 
position) calculated following iKovalev et alJ d2005h . The reso- 
lution limits were used as upper limits to the component sizes. 
Overall, nearly one-fourth of the cores were unresolved, namely 
98 at X-band and 87 at S-band. The projected linear sizes could 
also be calculated for the sources with known redshifts. The cor- 
responding distributions of the angular and projected linear sizes 
of the VLBI core components are presented in Fig. [TT] (middle 
and bottom panels), with median values of 0.28 mas (1.90 pc) 
and 1.04 mas (6.75 pc) at X-band and S-band measurements, 
respectively. 

The brightness temperature o f the core componen ts in the 
source rest frame is given by (e.g.. lKovalev et afll2005l) 



21n2S C01 - e J 2 (l +z) 

7lk $maj ^min 



(i) 



where k is the Boltzmann constant, S core is the VLBI core flux 
density, # ma j and 9 m i n are the FWHMs of the elliptical Gaussian 
component along the major and the minor axes respectively, A is 
the wavelength of observation, and z is the redshift. Measuring 
the flux density in Jy, component sizes in mas, and wavelength 
in cm, we obtained 



„ _ S core (l +z) 

l b, core - "Z ~ K-> 

"maj "min 



(2) 



where a A = 1.639 x 10 for X-band (A = 3.5 cm) and a A = 
2.337 x 10 11 for S-band (A = 13 cm). For sources without mea- 
sured redshifts (see Table [T) we used z — to determine a lower 
limit to the brightness temperature. Another fraction of sources 
that also yielded lower limits to the 7b values was a category of 
objects with unresolved cores. Altogether, about 40% of the T\, 
values are lower limits. 

The distributions of the brightness temperatures of the VLBI 
cores at 2.3 GHz and 8.6 GHz shown in Fig. [12] are broad cov- 
ering more than four orders of magnitude with the median value 
of 2.5 x 10 11 K at both frequencies, which is close to both the 
equipartition value of ~5 x 10 11 K (lReadheadlll994[) and the 
invers e Compton limit of ~10 12 K ( iKellermann & Paulinv-Tothl 
1 19691) . The higher brig htness temperatures may be attributed to 
Doppler boosting ( e.g.. IKovalev et al|2005t iHoman et al]l2006c 
ICohen et al.l 120071 IKellermann et al.ll2007h . The histograms in 
Fig. [121 obtained at different frequencies of 2 and 8 GHz, look 
similar because the maximum measurable 7b is determined by 
the physical length of t he effective projected baseline and th e 
core flux density (e.g.. IKovalev et al.l 120051 lKovalevll2009bl) . 
Therefore, for VLBI observations simultaneously made at dif- 
ferent frequencies, we would expect to obtain close 7b values, 
taking into account the flat radio spectra of the core components 
in most cases. The VLBI core parameters, such as the fitted flux 
density, sizes, or estimated resolution limits, and derived bright- 
ness temperature separately for 2.3 GHz and 8.6 GHz are listed 
in Table|7] The statistics of the distributions shown in Figs.l9UT2l 
are summarized in Table [8] The mean values of the distributions 
of 0maj, Trnaj, and 7\> were calculated u sing the ASURV survival 
analysis package dLavallev et al.llT992T) . 

3.5. Spectral properties 

Flux density and parsec-scale structure variations in AGN jets 
require the observations to be carried out simultaneously at dif- 
ferent frequencies in order to properly investigate the spectral 
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Source 


Epoch 


2.3 GHz 


8.6 GHz 






Score 






T b 


Score 




Smin 


T b 






Jy" 


mas 


mas 


K 


Jy" 


mas 


mas 


K 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


J0006-0623 


2002/03/06 


1.53 


1.87 


1.74 


1.48e+ll 


1.55 


0.81 


0.27 


1.56e+ll 


J001 1-2612 


1999/03/08 


0.40 


0.57 


0.57 


6.10e+ll 


0.33 


0.35 


<0.35 


>8.92e+10 


J0017+8135 


2002/03/06 


0.56 


1.00 


0.60 


9.60e+ll 


0.54 


0.41 


0.18 


5.16e+ll 


J00 19+7327 


1999/03/08 


0.68 


1.51 


0.84 


3.51e+ll 


0.14 


0.52 


0.52 


2.36e+10 


J0022+0608 


2002/01/16 


0.44 


1.19 


0.62 


>1.39e+ll 


0.41 


0.16 


0.16 


>2.55e+ll 


J0027+5958 


2001/10/29 


0.08 


0.82 


<0.82 


>2.62e+10 


0.08 


0.20 


<0.20 


>2.97e+10 


J0035+6130 


2002/12/11 


0.09 


1.93 


<1.55 


>6.97e+09 


0.09 


0.45 


0.45 


>7.29e+09 


J0050-0929 


2002/03/06 


0.56 


1.66 


0.92 


>8.61e+10 


0.74 


0.55 


0.19 


>1.19e+ll 


J0059+0006 


1999/05/10 


0.29 


1.41 


1.41 


5.83e+10 


0.10 


0.55 


<0.55 


>9.39e+09 


J0102+5824 


2002/01/16 


0.57 


0.96 


0.72 


3.17e+ll 


1.54 


0.09 


<0.07 


>6.89e+12 



Columns are as follows: (1) source name (J2000.0), (2) date of observations (YYYY/MM/DD), (3) fitted Gaussian core flux density, (4) FWHM 
major axis of the fitted core Gaussian, (5) FWHM minor axis of the fitted core Gaussian, (6) core brightness temperature, (7)— (10) the same as 
(3)-(6), respectively. Table |7]is published in its entirety in the electronic version of the Astronomy &■ Astrophysics. A portion is shown here for 
guidance regarding its form and content. 

Table 8. Statistics of distributions in Figs.lMHl 



Frequency, [GHz] 


S VLBI 


[Jy] 


S unre 


,[Jy] 


S unres/S'vLBI 


S core/ S VLBI 




Mean 


Median 


Mean 


Median 


Mean Median 


Mean Median 


2.3 


1.02 ±0.08 


0.76 


0.47 ± 0.03 


0.34 


0.51 ±0.02 0.51 


0.72 ± 0.02 0.75 


8.6 


1.04 ±0.10 


0.69 


0.54 ± 0.05 


0.32 


0.52 ±0.02 0.51 


0.73 ± 0.02 0.75 




S core i 


[Jy] 


^core j 


[mas] 


^core proj ; [pc] 


log(r b , core , [K]) 




Mean 


Median 


Mean 


Median 


Mean Median 


Mean Median 


2.3 


0.69 ± 0.05 


0.49 


1.18 ±0.09 


1.04 


7.56 ± 0.66 6.75 


11.68 ±0.05 11.40 


8.6 


0.76 ± 0.08 


0.47 


0.33 ± 0.03 


0.28 


2.18 ±0.24 1.90 


11.69 ±0.07 11.41 



properties of the sources. The RDV observations perfectly fulfill 
this requirement and enable us to study spectral characteristics 
of the observed sources by using (i) a detailed reconstruction and 
analysis of a spectral index distribution map of a source and (ii) 
applying statistical methods to derive the general spectral prop- 
erties and possible correlations/trends. In this paper, we imple- 
mented mostly the latter approach. For the statistical analysis, 
we used the calculated values of total VLBI flux densities from 
the 2.3 GHz and 8.6 GHz images (see TableO to derive the in- 
tegrated spectral index. The vast majority of sources have flat 
spectra, a > -0.5, for more than 95% of the sample (Fig. [TJl 
top). 

We also derived median spectral indices for the VLBI jets 
by applying the following procedure. First, for each source the 
X- and S-band images were obtained with the same pixel size 
(0.2 mas) and the same resolution by convolving them with the 
average restoring beam. A spectral index map was obtained by 
aligning the optically thin parts of the jet. Using the X-band data, 
the total intensity jet ridge line was then constructed with nearly 
equally spaced points separated by a distance of the pixel size. 
The procedure stops where the peak of a Gaussian curve fitted to 
the transverse jet profile becomes fainter than 4rms noise level 
of the image. The position of the core component taken from 
the model fit was chosen as a starting point. We note that the 
VLBI core is often but not always the brightest feature on the 
map. Thus, the image phase center, where the CLEAN algorithm 
tends to place the brightest component, would be an incorrect 
choice of starting point in such cases. Applying the ridge line to 
the spectral index distribution map, we extracted the values of a 



along the jet. After that, we masked out the core region within 
/"core < 0-5(^0 + d-p 5 , where b$ and are the beam and fitted 
core sizes, respectively, along a position angle </> of the inner jet. 
Finally, the jet spectral index was estimated as the median of the 
remaning values. The described method provided results for 319 
sources. The other 5 1 objects do not have enough structure, mak- 
ing the application of the procedure impossible. The measure- 
ments of the jet spectral indices are mostly distributed accross a 
range from -0.2 to -1.2 with the median value of a^_ g = -0.68, 
indicating optically thin radiation (Fig. [13] bottom). The spec- 
tral index of optically thin synchrotron radiation parametrizes 
the energy spectrum of relativistic radiative particles. Assuming 
a power-law energy distribution N(E) = NoE~ y , the power index 
y — 1 - 2a has the median value of y — 2.4 in the jets from the 
studied sample. The only positive value of a median jet spectral 
index was registered for the quasar J0927+3902, which has a 
short (~3 mas) jet with a resolved feature several times brighter 
than the core at both frequencies. This feature that dominates in 
the structure is optically thick at 2.3 G Hz and 8.6 GHz and h as a 
spectral index of 0.36. We also refer to lKovalev et al.l (1 19991) for 
the simultaneous total flux spectrum of the object. 

The spectral index on average decreases along the jet ridge 
line (Fig. [Pil l owing to the spectral aging caused by radiative 
losses of electrons. It slightly flattens in bright jet knots, which 
might be explained by shock acceleration at least in some cases. 
Using a linear least squares method, we estimated the spectral 
index gradients in the jets of 228 sources (Fig. IT5V This analy- 
sis was done only for sources with at least 10 ridge line points 
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Figure 13. Distributions of the 2.3-8.6 GHz spectral index of 
total VLBI scale emission {top), VLBI cores {middle), and jets 
{bottom). 
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Figure 15. Distribution of spectral index gradient along the total 
intensity ridge line for 228 sources. 



beyond the core region. The mean value of the spectral index 
gradient is — 0.06^2"2Z mas -1 , where the errors are given at the 
68% confidence level. 

Spectral indices for the VLBI cores were obtained as a value 
of the ff-map at pixels corresponding to the starting point of the 
ridge line (X-band core). The distribution of the core spectral 
indices (Fig. [13] middle) has the median value of orS°^ = 0.28, 
reflecting that the cores are partially opaque. 

We also examined the connection between total spectral in- 
dex and milliarcsecond compactness. The results for X-band 
data are shown in Fig. [16] the S-band dependence is qualitatively 
similar. The non-parametric Kendall's r-test confirms a positive 
correlation at a confidence level 99.9%. The explanation of this 
visible trend is that sources with steeper spectra are expected to 
have lower compactness indices, because a significant fraction 
of their VLBI flux should originate from optically thin extended 
jet emission. 
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Figure 16. Integral spectral index versus milliarcsecond com- 
pactness index at 8.6 GHz. The trend reflects the tendency for 
steep spectrum sources to be less compact. 
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Figure 17. 8.6 versus 2.3 GHz power-law indices of brightness 
temperature gradients with distance to the core {left) and with 
jet component size {right), dashed ellipse represents the lcr error 
areas of the k- and £-index distributions. 
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Figure 18. Distributions of a power index / assuming a d oc 
r dependence, where d is the transverse width of the jet at a 
distance r from its apparent origin at 2.3 GHz {left) and 8.6 GHz 
{right). 



3.6. Brightness temperature evolution along jets 

The sources with prominent jets were of particular interest to 
us, since they allow us to investigate the evolution of brightness 
temperature along the jet both as a function of distance to the 
core r, and the component's size d. Typically, a blob of rela- 
tivistic plasma, which is detected on a VLBI map as a bright jet 
component, loses a substantial amount of energy through syn- 
chrotron radiation and adiabatic expansion and quickly becomes 
dim propagating down the outflow and increasing in size. These 
factors lead to the rapid decrease in brightness temperature along 
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Figure 14. Spectral index distribution in J2202+4216 (left) calculated between 2.3 GHz and 8.6 GHz with the 8.6 GHz total intensity 
contours overlaid. White dots represent the total intensity ridge line, along which we plot the profiles of total intensity (top right) and 
spectral index (bottom right). The lowest contour is plotted at 0.18% of the peak brightness of 1395 mJy beam -1 . The shaded ellipse 
represents the FWHM of the restoring beam of 2.34 x 1.74 mas at PA = -14.0°, derived as the average between the corresponding 
interferometric restoring beams at 2.3 GHz and 8.6 GHz. The vertical orange dashed line shows the convolved radius of the VLBI 
core calculated along the inner jet direction. 



the jet. The evolution of 7t> can be well-fitted by a simple power- 
law functions 7b °c r~ k and 7b °c d~%. We analyzed the bright- 
ness temperature gradients for 30 sources with a rich jet structure 
consisting of at least three model fitted jet components at both 
frequencies. The power-law index k varies between 1.2 and 3.6, 
with the average value of ghz ~ ^2 GHz = 2.2 + 0. 1 . The power- 
law index ^ varies between 1 .4 and 4.4, with the average value of 
£s GHz = 2.7+0.1 and ^2 GHz = 2.6+0.1. Distributions of the fitted 
power indices k and £ are shown in Fig. [17] The power indices k 
and ^ can also be used to test jet mod els, e.g. adiabatic expansion 
in a shock-in-jet model suggested by iMarschej (I 1 990t) . 

Assuming that the width of a jet has a power-law dependence 
on distance as d oc r 1 , the power-law index I can be readily es- 
timated as k/%. The obtained distributions of I (Fig. fT8T > have 
median values of 1.2 at S-band and 0.8 at X-band, suggesting 
that jet regions probed by the RDV observations at 8.6 GHz 
are still collimating and therefore accelerating, while at 2.3 GHz 
the jets are more freely expanding and switch to a deceleration 
regime. Thi s is consistent with the positive parallel accelerations 
detected bv lPiner et al.l d2012l) analyzing te n-year jet kinematics 
based on the 8.6 GHz RDV data and also lHoman etalJ ( 120091) 
for 15 GHz MOJAVE data. 



3.7. Gamma-ray bright AGNs on parsec scales 

Among the 370 objects analyzed in this paper, 147 sources 
(40%, see Table [1} have been positionally associated with y- 
ray detections within their corresponding 95% confidence region 
made by the Large Area Telescope (LAT) on board the Fermi 
Observatory at en ergies above 100 M eV within two years of sci- 
entific operations (lAbdo et al.ll201 ll) . The list of the 147 sources 
include 99 of 251 quasars (40%), 39 of 46 BL Lacs (85%), 
6 of 31 radio galaxies (20%), and 3 of 42 optically uniden- 
tified sources (7%). A significantly higher detection rate for 



BL Lacs is attributed to their harder spectra, which allows them 
to be detected more easily by the Fermi LAT (lAckermann et al.l 
1201 ll) . We note that the vast majority of y-ray detections asso- 
ciated with AGNs have significant radio flux densities on par- 
sec scales as indicated by the highly successful rate of iden- 
tifications of the Fermi LAT objects with a r adio parsec-scale 
catalog of extragalactic jets (iKovalevi l2009al) . The Fermi era 
has already heralded a number of important links between y- 
ray emission and parsec-scale properties of AGNs. It has been 
shown t hat LAT-detected so urces are brighter and more lu- 
minous (jKovalev et al. 20091). have h igher apparent jet speeds 
(iLister et al.l l2009bl: iPiner et al.l 1201 2|) . higher Doppler factors 



(ISavolainen et aL 20 id), and characte rized by wider apparent 



opening angles ( Pushkarev et al.ll2009l) . There is also substantial 
evidence that y-ray emission likely originates at t he base of rel- 
ativis t ic jet close to parsec-scale radio core (e.g., Kovalev et alj 



120091: iJorstad & MarscheifcOlOtlPushkarev et al.ll2010l) 



We compared the flux densities of the VLBI core com- 
ponents for the LAT detected versus non-LAT detected RDV 
sources (Fig. IT9b and found that the y-ray bright sources have 
higher radio core fluxes, with a mean value of 0.89 Jy versus 
0.56 Jy at 2.3 GHz and 1.08 Jy versus 0.55 Jy at 8.6 GHz. In both 
cases, the Kolmogorov-Smirnov test indicates a probability p < 
0.001 that the corresponding samples are drawn from the same 
parent population. To complete a more detailed investigation of 
the y-ray/radio jet relation, we also examined the core brightness 
temperature 7b, core- Gehan's generalized Wilcoxon test from the 
ASURV survival analysis package indicates, at a > 99.9% con- 
fidence, that the 7b, core values for LAT-detected sources are sta- 
tistically higher than those for the rest of the sample at 8.6 GHz 
(Fig. [20] right), which is indicative of higher Doppler-boo sting 
factors. For the 2.3 GHz data, the difference is marginal (92.3% 
confidence) possibly due to the different spectral properties of 
the radio core and innermost jet components not being resolved 



14 



A. B. Pushkarev, Y. Y. Kovalev: Single-epoch VLBI imaging of active galactic nuclei at 2 GHz and 8 GHz 



2FGL 



2.3 GHz 



20 



□ 2FGL 



8.6 GHz 



20 



0.1 



0.1 

(Jy) 



Figure 19. Flux density distributions of VLBI core components 
at 2.3 GHz (left) and 8.6 GHz (right). The shaded areas represent 
LAT-detected (2FGL) sources. 
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Figure 20. Brightness temperature distributions of VLBI core 
components at 2.3 GHz (left) and 8.6 GHz (right). The shaded 
areas represent LAT-detected (2FGL) sources. 



10 ,„11 ,„12 ,„13 



with lower frequency o bservations. These fin dings confirm the 
early results reported bv lKovalev et al.l (12009b . 

Additionally, we examined the distribution of the median 
spectral indices for the VLBI jets obtained as described in § [ 



The distributions of aj et for 135 LAT-detected and 184 non- 
LAT-detected sources (Fig. l2"TT i differ at the significance level 
of 99.8% as indicated by a K-S test. The corresponding medians 
are clat.y = -0.60 and Q'lat.n = -0.72, which correspond to 
the power indices 7lat_y = 2.20 and 7lat_n = 2.44 of the en- 
ergy spectrum of relativistic radiative particles. This difference 
can be explained by a larger fraction of high-energy synchrotron- 
emitting electrons in y-ray bright AGNs being produced by ei- 
ther more energetic or more frequent flares. 



4. Conclusions 

We have obtained single-epoch images for a sample of 370 
AGNs based on the global VLBI observations carried out si- 
multaneously at 2.3 GHz and 8.6 GHz with the participation of 
ten VLBA antennas and up to ten additional geodetic radio tele- 
scopes. The parameters of the Gaussian model components rep- 
resenting the source structure have also been determined. The 
sample contains 251 quasars, 46 BL Lacertae objects, 31 ra- 
dio galaxies, and 42 optically unidentified sources. At least 97% 
of the sources manifest one-sided core-jet structures suggesting 
small angles to the line of sight. Even on milliarcsecond scales, 
the sources are highly compact and strongly core-dominated. 
Almost one-fourth of the core components are completely unre- 
solved on the longest baselines of the global VLBI observations. 

Of the 370 observed sources, 96% have flat-spectrum par- 
tially opaque cores with a median value of spectral index of 
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Figure 21. Distribution of the 2.3-8.6 GHz spectral index for 
the VLBI jets detected (top) and not detected (bottom) by the 
Fermi LAT. See discussion of the positive value found for the 
quasar J0927+3902 in Sec. [331 



ffcore ~ 0.3. The jet components are usually optically thin with 
a median value of c?j et ~ -0.7, which corresponds to y ~ 2.4 
assuming a power-law energy distribution N(E) - NqE~ 7 of ra- 
diative particles. The spectral index on average decreases along 
the jet ridge line due to the spectral aging. It slightly flattens in 
bright jet knots, which might be explained by shock acceleration 
at least in some cases. We have found that the estimated bright- 
ness temperatures of the VLBI cores are affected by Doppler 
boosting, their typical values are found at the level of 10 11 K, 
reaching up to 10 13 K in rare cases. We note that 40% of mea- 
surements yield only lower limits either because of an unknown 
redshift or unresolved core component. No systematic difference 
was registered for values of 7\, core obtained at 2.3 GHz and 
8.6 GHz. 

The evolution of the brightness temperature of the jet com- 
ponents as a function of distance to the core, r, and transverse 
size, d, was fitted with power-law dependencies of 7\, °c r~ k and 
7b oc d~^, respectively, for 30 sources with straight and promi- 
nent jets, for which the mean values of &8ghz ~ kiGHz = 2.2, 
^8GHz = 2.7, and ^2GHz = 2.6 can be used to test jet models. 
On the scales probed by the 8.6 GHz observations, a jet geom- 
etry for these sources shows evidence of further jet collimation, 
while on the larger scales probed by the 2.3 GHz observations, 
the outflows expand more freely, perhaps switching into the de- 
celeration regime. 

Among the 370 presented sources, 147 (40%) that are posi- 
tionally associated with high-confidence y-ray detections from 
the Fermi LAT Second Source Catalog have been found to have 
higher core flux densities and brightness temperatures suggest- 
ing preferentially higher Doppler-boosting factors. Furthermore, 
we have found that the LAT-detected AGNs are characterized 
by a less steep radio spectrum of the optically thin jet emission, 
most probably due to a larger fraction of high-energy radio emit- 
ting electrons being produced by flares. 

Our acquired VLBI images together with their respective 
Gaussian models may be used for future astrometric projects to 
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improve the accuracy of the source positions by taking into ac- 
count their milliarc second structure. The highly compact bright 
sources would be potentially the most suitab le targets for spa ce 
VLBI studies with the RadioAstron mission (lKardashevlll997l) . 
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